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By the application of the generalized Israel junction conditions we derive cosmological equations 
for the fourth-order f{R) brane gravity and study their cosmological solutions. We show that there 
exists a non-static solution which describes a four-dimensional de-Sitter (ciS'4) brane embedded in a 
five-dimensional anti-de-Sitter (AdSs) bulk for a vanishing Weyl tensor contribution. On the other 
hand, for the case of a non-vanishing Weyl tensor contribution, there exists a static brane solution 
only. We claim that in order to get some more general non-static f{R) brane configurations, one 
needs to admit a dynamical matter energy-momentum tensor in the bulk rather than just a bulk 
cosmological constant. 
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I. INTRODUCTION 



, Similarly as in the standard general relativity, it is interesting to consider generalizations of the brane universes 
Some of such generalizations are the higher-order brane gravity theories of which the simplest is f{R) gravity 
4 1 (for the most recent reviews see Refs. [y, Q). However, the combination of brane models with higher-order 
theories is non-trivial, since, except for the Lovelock (or, in the lowest-order, the Gauss-Bonnet) densities 0, 
one faces ambiguities of the quadratic delta function contributions to the field equations. This problem was first 
challenged succesfuUy in our earlier works d, in which we found the ways to avoid ambiguities not only for f{R) 
(— I ' brane theories (see e.g. [HIIi^Dj but also for more general actions which depend arbitrarily on the three curvature 
I— I, mvariants f{R,RabR°''',RabcdR°''"^'^), although the linear combination of these invariants was studied in Ref. One 
of the methods applied, was the reduction of the fourth-order brane gravity to the second-order theory by introducing 
an extra degree of freedom - the scalaron p, [I^l- Such a procedure leads to the second-order gravity which is just 
' the scalar-tensor Brans-Dicke gravity [13] with a Brans-Dicke parameter uj = 0, and an appropriate scalaron potential 
' (with the scalaron playing the role of the Brans-Dicke field) . We then obtained the Israel junction conditions [15j which 
T— I ■ generalized both the conditions obtained in Ref. [l^ [l3| for the Brans-Dicke field without a scalar field potential, 
' and also the conditions derived in Ref. 0, [l^ for f{R) brane gravity. The junction conditions which did not assume 
, scalaron continuity, but for a static brane, were presented in Ref. [20|. 

■ In this paper we apply the Israel junction conditions for /(i?) Friedmann- Robertson- Walker metric brane configu- 
rations and study the set of their admissible cosmological solutions. In Section II we present f{R) brane models and 

■ derive the set of field equations. In Section III we apply the field equations to cosmology. In Section IV we give our 
^ I conclusions. 

• ^ . 

^ , II. f{R) GRAVITY ON THE BRANE 

Let us consider the f{R) gravity on the brane described by the action Q 

1 

^5 JM, 



Sp = TT^ d^'xV^fiR) + Sbuik,p , (II.l) 

ZKr. 



where R is the Ricci scalar, Kg is a five-dimensional Einstein constant, Sbuik,p is the bulk matter action {p — 1,2), Mp 
is the spacetime volume. The action (jll.ip gives fourth-order field equations. It is then advisable to use an equivalent 
action 

Sp = [ d'xy^{f'{Q){R-Q) + f{Q)} (IL2) 
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where Q is an extra field (Lagrange multiplier), and f'{Q) = df{Q)/dQ. The equation of motion which comes from 
pi.2p is just Q = R, provided that f"{Q) ^ 0, so that f'{Q) may be interpreted as an extra scalar field (called the 
scalaron) 

= fiQ) = fiR) , (n.3) 

and the action can be rewritten as 

Sp = f d'xV^{^R-Vm + Sbuik.p- (II.4) 

where T^(0) = —(j)R{4>) + f{R{4>)) The action (jll.4[) is equivalent to a scalar-tensor Brans-Dicke gravity with 

a Brans-Dicke parameter lo — Q. One of the ways to derive the junction conditions for the theory described by the 
action (|II.4p is to append it with an appropriate Hawking-Lutrell boundary term which reads as [2ll | 

Shl, = -2(-l)^e / V^cl>Kd^x , (II.5) 

JdMp 

where K is the trace of the extrinsic curvature tensor Kab, h is the determinant of the induced metric hab — 9ab~^nanb, 
n" is a unit normal vector to a boundary dMp, and t — 1 {e = —1) for a timelike (a spacelike) brane, respectively. 
The total action of the theory is then 

Stotp = Sp + SLHp- (11.6) 

The variation of the action (III.6P leads to the following junction conditions for f{R) gravity theory (a, 6, . . . ~ 
0,1,2,3,5) 

W = , (II.7) 

S'^'^naHb - , (II.8) 

S'^'^hacUb = , (II.9) 

~hab[^^cn'']~[mhab + [4'Kab] = eKlS'^'Kahdb, (H.IO) 

where for an arbitrary quantity A we have defined a discontinuity (a jump) of A as: [A] = A'^ — A^ . Here Sab is the 
brane energy-momentum tensor (cf. later). In particular, the condition (jll.Tp comes from the requirement that the 
variation of the Hawking-Lutrell boundary term (|H.5P should vanish. 

These junction conditions can be compared with previously obtained in Ref. [l^ (see their Eqs. (12) and (13)) 
and in Ref. [11] (see their Eq. (3.11)). The difference is that we have not assumed the continuity of the scalaron on 
the brane. If we do so, i.e. assume that [0] = which due to the definition (|II.3p implies [i?] = 0, and additionally 
impose the mirror symmetry gab = 9ab{\n\), where n is a normal Gaussian coordinate originating at the brane, then 
the junction conditions (jII.7P - (III.10p take the form 

= , (H.ll) 

4>,cn' = -'^S'^'^hab , (H.12) 

Kab = ^{S^'^Kahdb - ^S^'Kd} ■ (11.13) 
Zip 4 

We can now express the equations above in terms of the Ricci scalar R instead of the scalaron (jj. Using the Gaussian 
coordinate system we obtain 

R = , (n.l4) 
K ^ , (H.15) 

f"{R)R.n = ~^S, (11.16) 

f{R)Kab = ^S^' , (H.17) 

where we have used the definition of a traceless part of the brane energy-momentum tensor 3°"^ = S'^'^hcahdb — 
{l/A)habS'^'''hcd, the definition of the traceless part of the extrinsic curvature tensor K""^ = K^'^hcahdb — / '^)habK'^'^hcd 
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and the definition of the trace of the brane energy-momentum tensor S = S°'^hab- The condition (jII.14[) is a conse- 
quence of the definition (|II.3p . The junction conditions (|II.14p - (|II.17p coincide with those obtained in Ref. [13] (Eq. 
(3.11)). 

Now we can apply the junction conditions (jll.TP - dll.lOp for f{R) gravity in order to obtain the effective Einstein 
equations on the brane. The manipulation of the Gauss-Codazzi equation ^22.] 

'^'^Rcdh^^n' = -D,K\ + DaK (11.18) 

where Da means a 4-dimensional covariant derivative on the brane, leads to a standard decomposition of a four- 
dimensional Einstein tensor ^■^'^Gab in the form 

= KKab - K^Kb. - \KbiK^ - i^-i^.) - ^^^Eab 

+ \[^''^G,dh\h\ + Kbi^^^^Cdn^n" - ^'^G)], (11.19) 

where Eab is an electric part of the bulk Weyl tensor projected onto the brane. We also assume that in the neighborhood 
of the brane, the normal vector field n°- fulfills the geodesic equations vf^V = (geodesic gauge). Using this last 
assumption, the following relations are derived (see the Appendix) 

(5)n0 = '^^^U(t> + Kn''Va(t>+{n''Vaf(t>, (11.20) 
h\h\V,^d^ - DaDb(t> + Kabn^Vc^. (11.21) 

Here Vc means the 5-dimensional (bulk) covariant derivative. Assuming that the matter in the bulk has the form 
of the 5-dimensional cosmological constant T^b = —gab^^^^, the variation of the action (III. 61) gives the following field 
equations in the bulk 

'"^'^-fc = -^9abV{^) + \gab'^''^n^~\^:ab+'^9ab^''^^. (11.22) 

Iq) <p (p ' (p 

(5)^ ^ _^XM. ^ _vt/(0). (11.23) 

0(p 

Substituting (III.lip - pi.13p . (|II.20p . (|II.2ip . and pi. 221) to pi.l9p . one obtains the efi'ective Einstein equations on the 
brane as 



where 



and 



l^DaDb^l) + ^''^'^^hab - ^^^Eab, (n.24) 
6 (p 2(p 



Qab = -^Xifab - ^fhab) + ^ffab - f'bTca + \Kb{f,df''' - ^f^) (11.25) 



fab Xhab + Sab- (11.26) 



This should be appended by the conservation law for the matter energy-momentum tensor on the brane (see Appendix 

DaS% = . (11.27) 

III. f{R) FRIEDMANN COSMOLOGY ON THE BRANE 

We restrict ourselves to the case of the matter in the bulk in the form of the cosmological constant. This allows to 
assume that the bulk spacetime is an Einstein space 

<'^C?ab = -''^Ae//5afc, (III.l) 
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where ^^^Ae// > is an effective cosmological constant. The five-dimensional hne element reads as 



ds-^ ^ -6 (^^, t)dr + a^{n, t) 



dn^ , (111.2) 



1 — kr^ 

where k — 0, ±1. The electric part of Weyl tensor Eab can be expressed in the following form [23[ 

^'^'>Eab = J'iUaUb + ^{hab + UaUb)]. (III.3) 

In the case with vanishing J^, we deal with a non-static Friedmann-Robertson- Walker brane (IIIL2|) embedded in an 
AdS^ bulk [iJ]. Moreover, the junction condition (|IL12p requires that the trace of the brane energy-momentum tensor 
vanishes S = T + AX = (the second bulk equation (IIL23|) sets the scalaron to be a constant because of the constancy 
of the curvature in the bulk). The assumption that the energy- momentum tensor of the matter on the brane is a 
perfect fluid 

Sab = (P + p)UaUb + pgab , (III-4) 

which fulfills the barotropic equation of state p = wp with the four- velocity vector u"" ^ 5q, gives the effective f{R) 
gravity Friedmann equations on the brane in the form 



2(j) 4(/) acj) 



(III.5) 



(p + p)(9p + p-8A), 



+ -2-^^ (III.6) 



1 f Ki 
3 



(p-l-p)(21p+13p-8A). 



Now the brane energy-momentum tensor conservation law (jII.27P 



^=-3(u; + l)- (III.7) 
p a 

can be integrated to give 

p = poa~'^^+'^ (111.8) 

The requirement that the trace of the brane energy-momentum tensor vanishes imposes a condition that the energy 
density of the matter on the brane is constant, i.e., 

4A 

^ = ^° = "(rT3^ = °°"'*- ("^-^^ 

Multiplying pil.6[) by three, and adding it to (jIII.5[) . we get one cosmological equation to solve (for simplicity we 
consider flat fc = models only) 

The Eq. (jll.23[) forces the scalaron (j> to be constant cj) = (pQ, as well. On the other hand, the eq. (jll.22[) leads to the 
following relation 

^-^.(^)A.,,. (III.11) 
200 00 " ^ ' 

The Eq. pil.lip shows that in the case of a constant scalaron, the term '^^^Ag// plays the role of a five-dimensional 
effective cosmological constant in the bulk. We can independently fix the value of the 0o by a choice of the shape of 
the function W(4i) near the brane using (III.23|) as 

00 = 1^-^ f-^^'^Aez/V (III.12) 
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Now for a fixed value of V(0o) = Vq, Eq. (jlll.lip gives 



Combining piI.lOl) with (IIII.lip . and assuming that w = — 1, we obtain 

-,2 



3(^ + -)=^'^Ae//. (111.14) 
The non-static solution of piI.14[) takes the form 

a = ao exp (Hot) , (III. 15) 

and it is consistent with the solution of (|III.5[) for = 0. For w ^ —1, the continuity equation pil.7[) implies a 
constant scale factor a{t) = ao, and the generalized Friedmann eqs. pll.5|) and pll.6|) become inconsistent. The 
solution pil.isp describes an embedding of a de-Sitter (dS'4) brane in an anti-de-Sitter {AdS^) bulk provided that 



(note that because of the assumption that ^^''Aeff > 0, we need Vq > 2k,'^'^^'>A). The Eq. (IIII.16P is a fine-tuning 
condition for the value of the bulk cosmological constant ^^^A and the potential V{(j)), which is responsible for the 
value of (po and Vq. The special case with 

Vo = 2kI'-^'^A (III.17) 

gives Hq ~ 0, and the solution pn.l5|) describes a static Minkowski brane which is a fiat analogue of the Einstein Static 
Universe. In fact, the condition piL17p is a special case of the fine-tuning relation piI.16|) . and can be interpreted as 
a necessary and a sufficient condition for the existence of a static brane in the model with w — ~1 and a vanishing 
Weyl tensor contribution = 0. 

If the Weyl tensor contribution is non- vanishing, i.e., if 7^ 0, it is then possible to embed a static Friedmann- 
Robertson- Walker brane pil.2p in a bulk with the cosmological constant '^^A only [^J]. In such a case, the solution 
of pil.ip for the metric pil.2p has the form ^ 

a'^in) = f{n), where f{n) = je^^^^'^l + 5e-^"°^'"\ , (III.18) 

b^n) = where e(n) = ^e^"«\''\ - Je^^^"!"! , (III.19) 

fin) 

with the brane at n = 0. Using the transformation of the metric components 

X /("-) , e^(n) e^(n) + S „„x 
/ " ^ ^ and -1^ ^ ; , III.20 

j + S f{n) f{n) [-y-dy 

which is equivalent to a rescaling of the coordinates, we obtain a Minkowski brane (for n — 0). We then compute 
the non- vanishing components of the electric part of the Weyl tensor and the corresponding term T at the brane 
(n = 0) as 

E\ = E\ = E\ = -E\ =F= ^-M_(5)A,,, . (III.21) 

For w 7^ -1, the solution of Eqs. (|nr7| . (UlLSj) and pIL6l) gives 

1 - 3M ± 2^/M MVq - 2kIK) 

9M^1 ' -^-""^^= sI^iA)^ >Q ("^-^^^ 

3 [kIX^ ' 



Note that pil.221) and pil.23p requires that (^0 > 0. The solution above describes a fiat static Minkowski brane which 
is a flat analogue of the Einstein Static Universe. In fact, from pil.2ip and pil.23p we have 

The condition pil.24p means that to support a fiat static Minkowski brane one needs to fine-tuned the values of the 
parameters 0o, Vo, '■^•'A, A, 7 and S. 
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IV. CONCLUSIONS 



In this paper we have studied brane universes within the framework of the fourth-order f{R) gravity theory. We 
apphed the junction conditions obtained in our earher papers d, [13] in order to get the set of the field equations 
which were apphcable to cosmology. We conclude that for the matter with a barotropic equation of state p = wp 
on the non-static Friedmann- Robertson- Walker brane (IIII.2[) embedded in a five-dimensional anti-de-Sitter [AdS^^) 
bulk (with vanishing Weyl tensor contribution = 0), and the matter in the bulk having the form of a cosmological 
constant '^^^A, there is only one case with w = — \ that possesses the solution in the form of the exponential evolution 
PII.ISP which is a four-dimensional de-Sitter {dSi) brane embedded in a five-dimensional anti-de-Sitter [AdS^) bulk. 
The case with the Friedmann-Robertson- Walker brane (jIII.2p embedded in a bulk with the cosmological constant 
^^^A and non-vanishing Weyl tensor contribution F ^ Q allows the solution in the form of the flat static Minkowski 
universe ()III.22|) only. The cosmological constant '^'^^A in the bulk implies the constant curvature and whence the 
vanishing of the trace of the brane energy-momentum tensor S. This is an extremely strong condition that makes the 
energy density constant. In conclusion, we claim that more non-static configurations on the brane are possible if we 
assume a dynamical matter energy-momentum tensor in the bulk. 
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Appendix A: Derivation of useful geometric formulas from Section |lTJ 



The formula pi.20p can be obtained as follows 




(A.l) 



where we have used 



nVVcVd^ = n^Vc{n''VdH) - n^(Ven^)(Vd0) 



(A.2) 




= h\h\V,{h^^Ve(^) + h\h\{V,n^)na{Vecl^) + 

+ h\h\iffndVcVe(j) = 

= DaDb(l) + Kadh\{n'Vc)(^ = DaDb<j3 + 

+ ifa6(n"Vc)<^. 



(A.3) 



It is also useful to prove that 



(A.4) 



Appendix B: Derivation of the brane energy-momentum tensor conservation 

The conservation law for the brane energy-momentum tensor can be obtained as follows. We take the covariant 
derivative of the left hand side of the equation pi.lBI) multiplied by and get 

Dam\) - K\Da<t> + ^DaK%. (B.l) 

Next, we use the Gauss-Coddazzi equation III. 181 together with the condition pi.Tp gives 

Dam\) - K\Da<j)-d^^'>)R,dh\n'. (B.2) 
Contracting (jll.22[) with the induced metric h\ and the normal vector one obtains 

^'''^R.dh^n'' = -\h\n''c^,,d ■ (B.3) 
<P 

After substitution of (|B.3p to (|B.2[) one gets 

Da{4>K%) = K%Da4> + h\n^4>,,a ■ (B.4) 

Applying (|A.4p to (|B.4|) we have 

Da{^K\) = K\Da4> + Db{n'V,4>)^K\Da^. (B.5) 
Substituting (|IlT2)) into (|R5|) . we obtain the relation 

Dam\) = -^D,S (B.6) 
16 

Now taking the covariant derivative of the right hand side of the equation (III.13P multiplied by (/>, we have 

Comparison of (IB.6|) with (jB.7p gives the conservation law of the brane energy-momentum tensor 

DaS\^Q , (B.8) 

as required. 
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